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ABSTRACT 


Adopting  two  physical  postulates  and  using  only  the  language  of 
eigenfunctions  and  eigenvalues,  we  express  the  underwater  sound 
field  of  a harmonic  point  source  as  a sum  over  orthonormal  discrete 
modes  plus  an  integral  over  a continuous  set  of  normal  modes. 
Restrictions  are  constant  water  depth,  no  absorption,  parameters 
varying  only  with  depth,  and  bottom  material  treated  as  a fluid. 

For  the  Pekeris  case,  orthogonality  of  the  continuous  modes  is 
demonstrated  (throughout  this  report,  orthonormality  is  with 
respect  to  a weighting  function,  the  material  density),  and  the 
integral  over  continuous  modes  is  proved  identical  to  the  branch 
line  integral  (BLl)  of  Ewing,  Jardetzky,  and  Press.  We  show  that 
farfield  contributions  to  BLI  sire  practically  negligible  except 
for  contributions  from  narrow  "spikes"  or  resonances  of  the 
integrand.  Conditions  for  existence  of  resonances  are  stated. 
Resonances  lead  to  virtual  modes  (discussed  by  Labianca  and  by 
Tindle  and  Guthrie)  that  decay  faster  with  range  than  a discrete 
mode  does.  Some  computer  results  by  Stickler  are  examined 
semianalytically  in  terms  of  discrete  and  virtual  modes;  agreement 
is  fairly  good.  We  find  that  BLI  also  has  a "pseudoresonance", 
yielding  a virtual  mode,  when  a discrete  mode  is  near  cutoff. 
Comments  on  Bartberger's  method  of  improper  modes  are  added. 
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INTRODUCTION 


We  consider  the  propagation  of  cw  single  frequency  sound  in  a 
medium  that  comprises  a water  layer  of  constant  depth  H and  a bottom 
half- space,  treated  as  a fluid.  The  sound  speed  c varies  only  with 
depth  z,  which  is  measured  vertically  downward  from  z=0  at  the  water 
surface.  The  material  density  p also  may  vary  wiAh  z,  but  at  most 
it  varies  via  finite  discontinuities  at  one  or  more  values  of  z.  An 
omnidirectional  point  source  of  angular  frequency  oj  is  at  depth  z=d 
on  the  z axis.  Horizontal  range  r is  measured  from  the  z axis. 

The  acoustic  velocity  potential  <p(  r,  z ) — a time  factor  exp(+ixot) 
having  been  suppressed— can  be  represented  as  a sum  over  a finite  set 
of  discrete  normal  modes  plus  another  contribution  expressible  as 
either  a branch  line  integral  (BLl)  or  an  integral  over  a continuous 
set  of  normal  modes.  If  H/\  (\  is  the  acoustic  wavelength)  Is  small 
enough,  the  discrete  set  is  empty.  If  certain  artificial  boundary 
conditions  are  given,  the  continuous  set  can  be  empty  but  then  the 
discrete  set  is  countably  infinite. 

In  the  sections  that  follow,  we  first  discuss  the  normal  modes  of 
the  medium  in  the  absence  of  any  source,  using  the  language  of 
eigenfunctions  and  eigenvalues.  Throughout  the  report,  effects  of 
absorption  are  neglected  as  are  all  other  dissipative  effects,  e.g., 
the  slight  acoustic  leakage  through  the  water  surface  to  the  air  above. 
We  "connect"  these  modes  to  the  point  source  described  above  and  show, 
for  the  special  case  of  a Pekeris  problem,  that  the  sum  of  discrete 
modes  is  identical  to  that  found  by  Pekeris,^  Ewing,  Jardetzky,  and 

p 

Press  (EJP),  and  others.  Continuing  with  the  Pekeris  problem,  we 

find  the  continuous  modes,  normalise  them,  prove  that  they  are 

orthogonal,  and  demonstrate  that  the  integral  over  all  continuous  modes 

2 

is  identical  to  the  BLI  found  by  EJP.  This  result  adds  support  to  an 


1 


argument  by  Stickler*  that  the  EJP  branch  line  integral  is  to  be 
preferred  over  that  of  Pekeris. 

Next,  the  EJP  integral  is  examined.  It  turns  out  that  for 
sufficiently  large  r,  i.e.,  in  the  farfield,  the  integral  practically 
vanishes,  unless  in  some  part  or  parts  of  the  range  of  integration  the 
integrand  displays  large  and  narrow  peaks  or  troughs.  The  effect  of 
such  a "spike"  is  much  the  same  as  that  of  an  extra  discrete  mode  being 

introduced,  which  suggests  that  Bartberger's  "improper  modes"^*^  can  be 
explained  thus.  Labianca  discussed  a closely  similar  problem0  and, 
following  quantum  mechanical  nomenclature,  referred  to  the  "extra 
modes"  as  virtual  modes.  We  use  his  terminology  and  test  the  effects 
of  such  modes  against  computer  results  published  by  Stickler.'*  Reasonable 
agreement  is  found. 

It  should  be  noted  that  just  after  this  manuscript  was  drafted  a 

7 

preprint  arrived  from  two  New  Zealand  scientists.  They  applied 
Labianca* s method  to  the  kind  of  problem  discussed  here  and — still  in 
accord  with  quantum  mechanics — designated  the  "spikes"  that  lead  to 
virtual  modes  as  resonances,  a term  that  we  also  adopt.  There  is 
considerable  but  far  from  complete  overlap  of  Ref.  7 and  our  analysis. 

I.  Two  Postulates  Concerning  Normal  Modes 

The  term  normal  mode  can  have  various  meanings . It  may  ?np»n  a 

pattern  of  motion  characteristic  of  a specified  medium  (which  is  our 

meaning  here)  or  a pattern  characteristic  of  a specified  medium  with 

a specified  source  in  it.  There  are  many  other  examples.  (See  a 

8 

discussion  by  Biot  and  Tolstoy.  ) We  might  use  the  term  "exotic  mode" 

(in  no  pejorative  sense)  to  mean  one  that  appears  to  be  a valid 
mathematical  solution  of  the  particular  problem.  Yet  at  first  glance 
and  often,  too,  at  second  glance,  the  exotic  mode  has  surprising 
behavior,  e.g.,  it  diverges  to  infinite  amplitude  as  z-*®.  We  claim 
no  proficiency  in  that  powerful  method  of  solution,  integration  in  the 


2 


Introduction,  the  mathematical  functions  c(z)  and  o(z 


represent  the  behavior  of  sound  speed  and  density  must  be 


reasonable.  Finite  discontinuities  sire  allowable  good  approxima- 


tions to  rapid  continuous  changes,  but  c(z)  and 


have  si 


field  is  nonnegligible . 


depth  z where  the  sound 


The  special  case  in  which  c(z)-*«  as  z-»*  may  be  open  to  debate 
(or  perhaps  it  is  covered  by  Ref.  9)>  but  we  note  that  the  resulting 
bottom  reflection  coefficient  appears  to  have  unit  magnitude  for  all 
angles  of  incidence.  A more  serious  example  of  misbehavior  is  the 
pseudolinear  positive  gradient  of  sound  speed:  c(z)=co(l-0(z-z 

fi>0,  z5zq.  At  certain  angles  of  incidence,  this  c(z)  leads  to  bottom 
reflection  coefficients  of  magnitude  exceeding  unity.  ® Some  of 
Bartberger’s  difficulties  with  exotic  modes  apparently  stem  from  the 
use  of  this  equation  in  depths  at  -aich  it  has  become  unphysical. 

Postulate  2.  Each  normal  mode  of  a medium  represents  a conceivable 


ttern  of  behavior.  The  set  of  ell  such  modes — discrete  plus 


continuous — must  be  normalizable,  complete,  and  closed  on  the 


OSzgco.  (The  additional  property  of  orthogonality  is  extremely 
convenient,  but  in  physical  problems  it  probably  exists  automatically 
or  is  achievable.) 


The  preceding  statements  -equine  that  each  mode,  separately,  must 
satisfy  the  wave  equation  and  all  pertinent  acoustic  boundary  conditions 
and  must  be  quadrat ically  integrable.  The  last  requirement  appears 
to  rule  out  Bartberger's  "improper"  or  "bad"  modes, which  diverge 
at  z=®.  However,  we  must  not  dodge  the  fact  that  inclusion  of  these 


modes  often  affects  the  total  solution  in  ways  confirmed  by  different 
methods  of  evaluation.  This  matter  is  discussed  later. 

Another  intent  of  this  postulate  is  to  rule  out  a priori 
imposition  of  additional  properties  on  any  single  mode,  e.g. , that  it 
must  be  a traveling  vave  or  a standing  wave.  Only  the  totality  of 
all  modes  determines  the  behavior  of  the  sound  field. 

H.  Eigenfunctions  and  Eigenvalues  of  Our  Problem 


For  a single  frequency  point  source,  the  acoustic  velocity 
potential  CD  is  governed  by  the  inhomogeneous  Helmholtz  equation 


V + x (z)?  = ~b (r-?o ) , 


,2,  ' at 
& (z)  = -5 

c2U) 


0 S z I ® 


The  right-hand  side  of  Eq.  1,  a three-dimensional  delta  function, 
corresponds  to  having  an  omnidirectional  point  source  at  z^.  Since 
we  assume  cylindrical  symmetry  about  the  z axis,  Eq.  1 reduces  to 


-l/d  \ (do\  3^0  ,2,  » 5(r) 


The  homogeneous  equation  corresponding  to  Eq.  3 is  separable.  We  take 
q>(r,z)=R(r)Z(z)  and  proceed  in  the  usual  way,  with  separation  constant 
Kp.  The  z equation  is 


~ + [k2(z i-K2]  Z(z)  = 0 
dz 
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Each  acceptable  solution  Z nust  be  normal  1 sable  over  0£z£®  and  must 
obey  all  pertinent  acoustic  boundary  conditions  at  an  finite  depths. 
These  conditions  are  fl ) continuity  of  acoustic  pressure  and  therefore 
of  pZ  and  (2)  continuity  of  normal  particle  velocity  and  therefore 
of  dZ/dz.  At  the  water  surface  the  sole  condition  is  that  Z=0.  Apart 
from  the  requirement  of  nornalizability,  we  postpone  any  comment 
about  the  radiation  condition  at  infinity.  A pertinent  fact  of 
geophysics  is  that  c(z)  and  p(z)  will  always  exceed  their  maxima 
values  in  the  water  layer,  at  and  below  a sufficiently  great  depth  in 
the  bottom. 

It  follows  that  Ea.  1 nust  be  solved  as  an  eigenvalue  problem. 

The  set  of  acceptable  Z's  includes  a subset  (possibly  empty)  with  a 

2 

finite  number  of  discrete  Z 's,  marked  by  discrete  values  k of  the 

2 n n 

separation  constant  \ . The  rest  of  the  set  is  made  up  of  the 

continuous  Z's,  in  which  adjacent  Z’s  and  their  values  of  K2  differ 

infinitesimally.  For  them,  we  write  K2  as  k2. 


As  3iot  and  Tolstoy  have  shown,  orthonormality  for  any 
acceptable  Z's  nust  be  with  respect  to  a weighting  function  p(z). 
That  is,  if  K and  K*  correspond  to  two  values  of  K2,  th=  integral  to 
be  treated  is 


C(K,K'  ) = f p(s)Zx(z)ZK*(z)dz  = 6k  k,  or  5(K-K*  ) 


Following  postulate  2 of  Section  I,  we  Sosume  that  the  whole  set 
(Z(k,z))  is  complete  (in  practical  terms  this  assumption  means  that 
we  have  taken  pains  to  find  all  acceptable  Z's)  and  is  also  closed. 
Completeness  allows  writing  <?(r,z)  as 


[r,z)  = S.Ji(K,r)z(K,r) 


iiiaattKMiiWiwtwMi, 


Here  O implies,  in  effect,  a Stieltjes  integral — a sum  over  discrete 
Zn's  and  an  integral  over  continuous  Z's.  Closure  of  the  set  is 
expressed  by 

SKZ(K,z)z{K,zt)  = 8(z-z* ) . ( 

In  Eqs.  6 and  7,  it  is  assumed  that  the  Z's  have  already  been 
normalized.  Taking  z'=d,  the  source  depth  and  substituting  Eqs.  6 
and  7 in  Eq.  } leaa  to 


It  follows  that  (...)  vanishes  identically.  The  resulting  solution  is 

A(K,r)  = (£)z(K,d)H£2)(Kr)  , {[ 

in  which  the  proper  Hankel  function  has  been  chosen  to  ensure  an 
outgoing  cylindrical  wave.  The  final  expression  for  tp  is 

q>(r,z)  = iC  S^Z(K,d)z(K,z)H^2^(Kr)  . (l) 

The  arbitrary  constant  C,  allowed  because  this  problem  is  linear,  has 
subsumed  the  factor  l/k  of  Eq.  9*  The  Z's  in  Eqs.  7 through  10  are 
normalized . 

Orthogonality  is  readily  demonstrated  for  the  discrete  modes, 
provided  that  p(z)  is  piecewise  constant,  (if  this  behavior  is 
violated,  the  acoustic  wave  equation  must  be  reexamined;  the  compli- 
cations resemble  those  in  nonlinear  acoustics.)  We  treat  the  case  of 

a single  discontinuity  in  p(z),  at  z=H;  generalizing  to  any  finite 

2 

number  of  discontinuities  is  simple.  Writing  Eq.  k for  Z and  K , we 


multiply  through  by  Z^,  repeat  th17  process  with  n and  m interchanged, 
subtract  one  equation  from  the  o^her,  and  obtain  (with  primes  for 
d/dz) 


Z Z"  - Z z"  = ■T-  (z  Z'-Z  Z')  = (k  2-k  2)z  Z 
mn  nm  dz'mnnnr  ' n m ' m n 


After  multiplying  through  by  p(z),  we  integrate  both  sides,  from  z=0 
to  H and  from  H to  <»: 

/CO 

p(z)Z  Z dz  . (12) 

v ‘ m n ' ‘ 


The  Z's  vanish  at  the  pressure- release  surface,  z=0.  They  also 
vanish  at  z=°°  if  postulate  2 and  the  comments  about  c(z)  at  great 
depth  (see  the  paragraph  following  Eq.  4)  are  accepted.  The  left-hand 
side  of  Eq.  12  becomes 


‘«„Zmw<K>Z’ „b<H»  ‘ 


Here  subscripts  w and  b mean  "in  the  water,"  "in  the  bottom,"  and  the 

second  [ ] is  the  same  as  the  first,  with  n and  m interchanged.  The 

acoustic  conditions,  or  the  continuity  of  oZ  and  of  Z'  at  H,  ensure 

that  all  of  expression  13  vanishes,  and  hence  Z and  Z are  orthogonal 
2 2 m n 

for  kn  /k  . The  proof  just  given  does  not  imply  that  orthogonality 


must  fail  when  p(z)  varies  continuously. 


Still  considering  only  discrete  modes,  we  apply  Eq.  10  to  the 

Pekeris  problem.^  c(z)  and  p(z)  are  piecewise  constant,  being  cv>Py 

in  the  water  and  c,  >c  , n>p  in  the  bottom.  The  acceDtable  unnormalized 
i)  v’  v „ j 

Z's  are  evident  from  Eq.  4 (with  K=k  ) and  postulate  2: 


7 


..I  ..Cl Mm  Hurt.  Iuh  Cl  hilt . « J iW«  K. 


Bf  1 

P { 


Considering  only  CKd<H,  we  substitute  Eqs.  14  and  18  in  Eq.  10, 

and  then  repeat  with  Eqs.  15,  17,  and  18.  The  two  results  are 

identical  respectively  to  Eqs.  A48  and  A50  of  Ref.  1,  if  we  take 

C=-«p  . They  are  also  identical  to  the  corresponding  equations  of 
2W 

EJP.  There  is  nothing  novel  about  this,  for  no  one  has  ever  doubted 
such  agreement;  our  main  point  in  this  exercise  is  to  establish  the 
proper  value  of  C for  later  use.  There  is,  however,  a touch  of 
historical  irony.  Temporarily  adopting  the  language  of  eigenfunctions 
(see  Ref.  1,  pp.  53-55),  Pekeris  expressed  concern  over  the  failure  of 
these  functions  to  be  orthogonal  unless  PW=P^-  The  cause  was  overlooking 
the  need  of  a weighting  function  (Eq.  5).  The  irony  was  that  his  method 
of  integration  in  the  complex  plane  had  automatically  produced  the 
effect  of  proper  orthonormal  eigenfunctions  in  the  sense  of  Eq.  5. 


We  observe  that  Eq.  l6  allows  no  discrete  mode  unless  > it/2  and 
that,  in  conjunction  with  Eq.  15,  it  allows  at  most  a finite  lumber  of 
modes . 

III.  Continuous  todes  of  the  Pekeris  Problem 

Through  the  rest  of  this  discussion  we  consider  only  the  Pekeris 

problem,  without  absorption.  For  the  continuous  modes,  Eq.  4 still 

holds.  We  now  write  Z(k,z)  and  K=k  instead  of  kR,  with  k a continuous 

variable.  Equation  14  also  holds  with  these  changes  of  notation,  and 

the  expression  for  (3,  is  unchanged.  However,  Eq.  15  is  different 
2 2 X 

because  now  k <k^  (see  the  remarks  immediately  following  Eq.  15). 
Acceptable  unnormalized  Z's  in  the  bottom  are 
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with  n and  v being  two  functions  of  k yet  to  be  determined.  (The 
only  other  conceivable  form  for  Z would  be  a complex  exponential,  but 
then  the  acoustic  conditions  at  z=H  could  not  be  met.)  After  writing 
those  conditions,  we  obtain  two  equations  from  them.  One  is  the 
eigenvalue  equation, 


(^l)  tan^lH^  = (§)  tan(p2H+v)  > 


which — thanks  to  the  leeway  afforded  by  the  extra  constant  v — can  be 

2 2 ? 
satisfied  for  any  k <k^  . The  other  equation  yields  4 : 


■ (e;)  '“’V) + (^)  sl“2( V)  • 


We  normalize  Z(k,z)  to  obtain  IT^Z,  as  in  Eq.  5,  by  the  method  of 
"box  normalization."  That  is,  the  integral  from  H to  “ is  carried  out 
from  H to  (H+L)»H,  and  eventually  L is  allowed  to  increase 
indefinitely.  The  result  is 


rH  2 2 rE+L  2 

= J sin  (frjz)dz  + p^u' J sin  (fr,z+v)dz 


= const,  x L + 


(set)  "-I  si“(2vt2v)|i 


• (25) 


The  sine  term  is  also  const,  x L , v drops  out,  and  the  outcome  is 


To  demonstrate  orthogonality  in  the  sense  of  Eq.  5,  we  let  p^, 

32,v  correspond  to  one  value  of  k and  to  another  value, 

k' . (This  Pg  not  the  same  one  that  is  used  in  Eq.  15.)  Temporarily, 
we  drop  subscript  2 from  since  P^  does  not  appear  anywhere  in 

this  calculation.  Also  we  express  p and  p*  as 


« mn 

p = r 


p*  = 


(25) 


where 


m / m* 


m,m*  = 0,1,2, 


(26) 


This  step  restricts  P and  p'  to  discrete  changes.  However,  the  smallest 
change  is  it/L,  which  approaches  a differential  quantity  when  L is 
allowed  to  increase  without  limit. 


We  can  now  make  use  of  the  method  in  Eqs.  11  through  1?  for 
arbitrarily  large  but  finite  L because  the  continuous  modes  have  been 
rendered  discrete  by  Eqs.  25  and  26.  The  only  change  is  that,  in 
Eq.  12,  the  upper  limit  « is  replaced  by  H+L  on  both  sides  of  the 
equation.  Consequently,  sane  extra  terms  appear  in  the  left-hand  side 
of  Eq.  12: 


pJZjH+L)  Z;(H+L)  - Zq(H+L)  Z;(H+L)j  = 

X |p*  sin[p(H+L)+v]  cos[P’ (H+L)+v' ] - L CoslP(H+L)+v]  .sin[P'(H+L)+v* ] J 
1 127) 

The  right-hand  side  of  Eq.  27  results  from  Eq.  19.  With  the  use  of 
Eq.  2 k,  the  term  N^. . .p^  in  Eq.  27  becomes  simply  2/L.  We  remove  PL 
and  P'L  from  within  { } by  invoking  Eq.  25  and  noting  that 

sin(nm+0)  = (-)m  sin  9 cos(mn+0)  = (-)m  cos  0 . (23) 

The  outcome  is  that  the  right-hand  side  of  Eq.  27 — the  extra  terms  in 
Eq.  12— is 
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vc1  s u puli'  Wj.  ■ j&tk Sift 


1rr^Ti*:faTi 


x [0*  sin(0H+v)  cos(0'H+v' ) - 0 cos(0H+v)  sin(0'H+v')J  . (29) 


For  any  given  pair  of  modes  denoted  by  0 and  0',  this  expression 
vanishes  as  L-»».  Simultaneously,  the  upper  limit  H+L  on  the  integral 
in  Eq.  12  once  more  becomes  ®.  Therefore  the  continuous  modes  are 
mutually  orthogonal.  The  orthogonality  of  any  discrete  and  any 
continuous  mode  can  be  proved  similarly. 

If  postulates  1 and  2 are  observed,  finding  the  normal  modes  for 
more  complicated  cases  would  differ  mainly  in  detail.  Equation  10 
remains  correct,  and  for  the  discrete  modes,  so  do  the  steps  in 
Eqs.  11  through  13.  Normalizing  these  modes  would  be  accomplished  by 
direct  integration.  The  corresponding  treatment  of  continuous  modeB 
could  proceed  via  box  normalization.  However,  the  calculations  will 
be  more  difficult  than  they  are  for  a Pekeris  case. 

IV.  Equation  10  and  the  Branch  Line  Integral 


The  continuous-mode  part,  cpc>  of  <p  in  Eq.  10,  can  now  be  written 
as  follows: 


Pc(r.z)  = ic£  Nk2Z{k,d)Z(k,z)H^(kr)  , 


with  the  unnormalized  Z*s  fro*.  Eqs.  l^f  and  19  as  pertinent.  (Thanks 
to  Eq.  25,  and  all  functions  of  them  are  determined  by  the 

integer  m;  hence,  summation  is  correct).  As  L-»“,  the  changes  5n 
quantities  corresponding  to  Am=  ±1  are  essentially  infinitesimal  and 
the  sum  can  b°  replaced  by  an  integral  ever  m.  The  use  of  Eqs.  19  and 
25  shows  that 


, - Hi'iirt! , it*  li 


^ (i 


We  take  C=-«p  , as  was  found  for  the  discrete  Pekeris  modes  in 
w’ 

Section  II,  employ  Eqs . 21  and  2k,  and  for  the  case  0<(d,z)<H, 
obtain 


cp  fr,z)  = 2i/~ 


pw\  r P2  sin(^id)  sin(pxz)  H^2'(kr)  kdk 
^7  ^ COS2(P1H)  + (Py^/Pb^  Sin2^lR) 


It  remains  only  to  find  the  limits  on  k or,  more  precisely,  the  proper 

contour.  Equation  19  shows  that  k starts  at  k^  and  decreases  along 

the  real  axis  until  k=0.  Nothing  prevents  k2  fr^m  -chen  becoming 

negative  because  the  separation  constant  K2  in  Eq.  k could  have  been 

taken  with  either  sign.  Hence,  integration  continues  along  the 

(2) 

imaginary  axis  of  k.  The  negative  axis  must  be  used  because  '(kr) — 

originally  chosen  in  Eq.  9 to  ensure  outgoing  cylindrical  waves  as 
r-*»  —diverges  as  +i  |k|->«  . The  contour  of  integration  in  Eq.  32 

is,  therefore,  k.  to  0 to  -i«.  Equation  32  is  then  identical  to  the 

D 2 3 

branch  line  integral  found  by  EJP  and  preferred  by  Stickler  over  the 

corresponding  integral  of  Pekeris.1 

For  the  case  d<H  but  z>H,  sin(P^z)  of  Eq.  32  is  replaced  by 
P sin(pgZ+v): 


q>  (r,z>H)  = 2i( 


/ Pw\  ^ sin^d)  sin(pgz- 

\V  \ [p/  cos201H)  + (pw 


,z+v)  H^2'(kr)  kdk 


P2/Pb)2  sin2(  PjH)  j1/2  ' 


which  is  identical  to  the  result  of  EJP.  Because  there  is  no  longer 
any  distinction  between  "branch  line  integral"  and  "integral  over 
continuous  modes,"  we  use  the  abbreviation  BLI  for  the  right-hand  sides 
of  Eqs.  32  and  33. 
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V.  Some  Properties  of  the  Branch  Line  Integral 

Over  the  interval  OS  k<k^,  the  functions  f(k)  multiplying  H^2^ 
in  Eqs.  32  and  33  are  real  and  continuous  with  continuous  derivatives. 
Since  kdk  becomes  -|k|d|k|  on  05kS-i»,  similar  behavior  occurs.  As 
k-»k^ on  the  real  axis,  ft,  goes  to  zero  and  ordinarily  f(k)  vanishes. 
However,  in  the  special  case  that  P^H  is  an  odd  multiple  of  «/2  at 
k=k^,  f(k)  must  be  tested  for  a possible  singularity.  It  turns  out 
that  f(k)  remains  finite  at  k^  (see  Appendix  A).  This  case  corresponds 
to  having  the  highest  discrete  mode  occur  exactly  at  cutoff.  Both 
integrands  vanish  at  k=0  because  H^2^  has  only  a logarithmic 
singularity  there. 

For  any  real  k>0,  the  asymptotic  form  of  H^2^(kr)  at  large  r is 
a cylindrical^  spreading  and  outgoing  wave . For  -ilk  1^0,  H^2^  falls 
off  asymptotically  as  (|k|r)  ' exp(-|k|r).  7 o consider  the  behavior 
at  large  z,  we  must  use  Eq.  33.  For  any  Pg/O,  i.e.,  kC}^,  the  quantity 
sin(P2z+v),  regarded  as  a function  of  k,  oscillates  ever  more  rapidly 
as  z is  made  larger;  therefore  the  integral  approaches  zero  as  z-»®. 
Hence,  the  radiation  condition  is  satisfied  for  large  r and  z,  which  may 
have  seemed  doubtful  because  only  standing  waves  in  z,  not  downward 
progressing  ones,  enter  Eq.  19. 


Eq.  33  (z>H)  will  not  be  considered  further,  but  Eq.  32  is 
needed  to  make  a comparison  with  seme  of  Scickler's  results.  For 
sufficiently  small  |k|  on  either  axis  the  integrand  in  Eq.  32  usually 
approaches  M|k|H^2^,  with  M being  a constant;  the  first  correction 
terms  sire  in  |kp.  (For  an  exception  see  Note  at  the  end  of 
Appendix  B.)  At  and  near  this  limit  the  integral  is  known  in  exact 
form.'*"1'  When  k is  real  and  r exceeds  zero,  the  result  is 


(2 


*•  - %)\7)L  (kr)  ° (kr>  d(kr) 


On  the  negative  imaginary  axis,  with  k=-i|k|  and  r>0,  we  substitute 
5=-i|k|  and  find 


<p_ 


I)  (-iW42)<-‘Wr)  - (i)j 


( 


Therefore  the  quantities  -(2i/*r)  cancel  when  Eqs.  34  and  35  are  added 
to  integrate  "around,  the  corner”  near  k=0.  The  surviving  terms  in  ^ 
vanish  as  k-*0,  by  mutual  cancellation. 

If  in  Eq.  35  the  linear  approximation  M|k|  holds  for  |k|  large 
enough  to  ensure  that  |k|r£2jc,  then  Ih^2^-! |k|r) | is  not  larger 
than  about  6x10  and  decreases  rapidly  for  larger  |k|r.  In  contrast, 
at  real  kr»2«,  the  corresponding  tern  |H^2'(kr)|  is  about  0.3.  Apart 
from  its  role  in  canceling  the  term  (2i/jtr),  the  integral  along  the 
negative  imaginary  axis  is  obviously  negligible  for  |k|r>2jt.  The 
behavior  of  BLI  in  Eq.  32  then  rests  primarily  upon  the  integral  from 
kj^  to  k=0. 


It  is  most  unlikely  that  the  integral  in  Eq.  32  could  be  evaluated 
without  numerical  procedures.  Still,  additional  information  about  its 
properties  can  be  gained  from  analytical  approximations.  Obvious  first 
comments  are  that  the  factor  k in  f(k)  implies  only  small  contributions 
near  k=0  and  that  phase  oscillations  of  the  Hankel  function  will  tend 
to  limit  ;he  complete  result. 

Taking  the  negative  of  this  integral  to  compensate  for  inverting 
the  limits  to  0 and  k^,  we  divide  the  real  k axis  at  points 


Here  m^  and  a^,  chosen  anew  for  each  interval,  are  always  labeled  by 
the  higher  value  of  l.  This  polygonal  approximation  for  f(k)  is 
made  continuous  across  each  k^  by  requiring  that 

f(k t)  “ Vi  + a|  = mi+iki  + al+ 1 * (57) 

Ignoring  the  multiplying  factor  -2i(pw/pb)  that  precedes  the  integral, 
we  have 


f(k)H^2  V)dr  = £ jP  = E(J?+J?  ) > 

*/o  til ' "I  ;i  ' 


(kr)H^2)(kr)d(kr) 


42)(kr)d(kr) 


and 


(38) 

(39) 


(to) 


has  the  same  form  as  the  integral  in  Eq.  34;  the  whole  tern  is, 
therefore, 


I 

I 


s 


Jk^^r)  - kf_142)(Vir)]  * (4l) 
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There  is  an  exception  when  k^  ^=kQ=0,  for  then  the  last  term  in  Eq.  4l 
is  replaced  by  -(2i/nr)  as  in  Eq.  35.  However,  this  quantity  is 
canceled  by  a like  contribution  from  the  integral^  over  imaginary  k 
(Eqs.  3^  55)  and  will  be  ignored.  Handling  and  combining  it  with 

is  more  complicated;  we  relegate  the  details  to  Appendix  B. 


The  outcome  (Eq.  B13)  is 


V1 


- * kai/2 (mi+i"mi)kii/2  •’1V 

(42) 


i= l 


Evidently  the  main  contributions  of  BLI  to  the  mode  sum  will  stem  frcm 

f(k)  having  a slope  of  large  magnitude  near  k=kfe  and/or  from  parts  of 

the  finite  sum  being  large.  The  term  m.  ,-m.  for  the  polygonal  approxi- 

* x 2 

mation  to  f(k)  is  roughly  proportional  to  d f/dk  . Thus,  large 

contributions  to  the  mode  integral  cane  from  intervals  £k  in  which  the 

slope  of  f is  changing  rapidly.  For  example,  suppose  that  in  some  part 

of  0<k<kb  the  function  f(k)  has  a high  narrow  "spike"  which  is  almost 

a triangle.  The  polygonal  approximation  might  then  employ  only  two 

intervals:  M to  f and  I to  1+1  j m,  would  be  large  and  positive  and 

1 -1/2  , , 

m^+^,  large  and  negative.  Over  these  two  intervals,  k^  ' would  change 

only  slightly.  If,  in  addition,  r were  not  too  large  so  that  exp(-ik^r) 

did  not  change  greatly,  this  contribution  to  BLI  would  be  a virtual 

2 

mode  with  an  eigenvalue  kf  . As  r increases,  however,  both  the 
factor  r~^/2  and  the  phase  oscillating  exponential  term  would  make  the 
virtual  mode  become  rapidly  less  important. 

VI.  Locating  Resonances  and  Virtual  Modes 

As  we  have  said  before,  eventually  BLI  must  be  evaluated  numerically. 
However,  it  is  simple  to  name  approximate  locations  of  the  resonances 
described  in  Section  V.  From  Eqs.  14,  19,  and  32,  we  have 


k02  sin(p..d)  sin(&z) 

f(k)  = —5 5 P p , and 

Px  cos  (&jH)  + (Pwe2/Pb)  sin  (PjH) 


• M 


Only  real  k,  O^kS  k^,  is  considered.  The  constant  ( Pw/ ) in  the 
denominator  would  usually  be  some  0.25  to  0.5.  Because  Eq.  52  holds 
only  for  (d,z)<  H,  p^d  and  p^z  in  the  numerator  will  be  smaller  than 
PjH.  "lie  product  of  the  two  sine  functions  will  oscillate  (unless  d=z). 

For  k far  enough  from  zero,  the  product  can  change  sign;  on  the  average, . 
it  does  little  else.  For  sufficiently  small  k,  P^  and  pg  are  approximately 
k^  and  k^,  respectively,  and  the  denominator  remains  nearly  constant.  The 
magnitude  of  the  numerator  initially  rises  roughly  in  proportion  to  k; 
for  larger  k,  it  may  develop  some  oscillations  about  this  linear  rise. 

As  k increases  enough  so  that  P^  and  Pg  are  no  longer  nearly 

constant,  the  denominator  begins  to  govern  f(k).  For  sufficiently  large 

k^H,  the  angle  P^  goes  through  whole  cycles  in  which  the  cosine  (sine) 

term  vanishes  when  P^H  is  an  odd  (even)  multiple  of  it/2.  For  seme  range 

of  intermediate  k,  nothing  exciting  happens  because  the  oscillations  of 
2 2 

cop  and  sin  roughly  balance  each  other.  As  k comes  close  enough  to 

k.  , however,  p,  and  especially  P_  change  more  rapidly  with  k.  P-.  can 
^ 2 2 l/2  ^ 

never  be  less  than  (ky  -k^  ) ’>0,  but  Pg  approaches  zero.  In  this 

range  of  larger  k,  when  P1H=^2n+l)(n/2),  the  cosine  tern  becomes  very 

small.  The  sine  term  itself  is  nearly  unity  but  it  is  multiplied  by 

2 

a tern  smaller  than  Pg  . Hence,  the  whole  denominator  becomes  small 
and  j f ( k ) [ becomes  large  rather  abruptly.  As  k and  p^H  pass  through 
this  critical  region,  | f (k) | decreases  again,  about  equally  abruptly. 


It  fo.llows  that  the  locations  of  large  resonances  in  f(k),  with 
consequent  virtual  modes,  are  determined  by  two  simultaneous  conditions: 
(l)  k being  not  too  far  from  the  branch  point  k=k^  and 
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(2)  (3^H=(2n+l)it/2.  Moreover,  the  resonance  vill  be  stronger  and 

narrower  the  nearer  condition  2 gets  to  k=k^.  Conditions  1 and  2 

are  necessary  but  not  sufficient  for  appearance  of  a resonance 

because  the  magnitude  of  the  numerator  night  be  very  snail  where 

resonance  would  otherwise  have  been  predicted.  We  c.ote  that 

substantial  resonances  cannot  occur  on  the  imaginary  axis  of  k because 

remains  large  everywhere.  It  is  interesting  to  observe  that 

conditions  1 and  2 are  also  the  requirements  for  explaining  the 

12 

bottomborne  first  arrival. 

Even  for  fairly  low  frequency  ' ^ shallow  water,  a few  resonances 
are  possible.  As  frequency  and  water  depth  increase,  i.e.,  as  H /k 
increases,  the  number  of  potential  resonances  rises.  An  approximate 
upper  bound  is  the  number  of  integers  in  the  range  k H/n  to 

^w2-\2)1/2h/«- 

VII.  Application  to  Stickler’s  Figure  11 

Stickler  has  published  numerical  calculations  of  transmission 

loss  versus  range  for  different  sets  of  assumed  conditions.^  His 

Fig.  8 approximates  a Bekeris  case  in  shallow  water,  with  H=150  ft, 

d=20  ft,  and  z=40  ft.  There  is  a slight  decrease  of  c(z)  with  depth 

in  the  water  layer  but  there  is  a constant  and  a density  ratio  of 

1.4.  At  a frequency  of  100  Hz,  there  is  one  discrete  mode.  Plotted 

against  range  r from  less  than  1000  ft  cut  to  21,000  ft.  Stickler's 

Fig.  11  shows  transmission  loss  at  100  Hz.  Tie  figure  gives  separately 

the  contribution  of  the  one  discrete  code  and  tie  contribution  of  the 

discrete  mode  plus  BLI.  Moreover,  as  a matter  of  particular  interest 

for  the  present  analysis,  Bartberger  had  contributed  20  points  spaced 

along  the  discrete- plus- continuous  curve  and  computed  by  his  method  of 
if*  5 

improper  modes.  ’ His  agreement  with  Stickler's  plotted  curve  is 
eveiywhere  good. 


Since  postulate  2 (Section  i)  rules  out  improper  modes,  we  seek 
another  explanation  of  this  agreement  in  terms  of  Sections  V and  VI. 
The  first  step,  approximating  Stickler’s  Fig.  8 by  a Pekeris  problem, 
required  the  estimating  of  a constant  average  c^  in  the  water  layer. 
Otherwise,  his  data  were  retained.  At  100  Hz,  the  chosen  values  were 

v 

k^=0. 125916  and  k^= 0.12 4223.  The  eigenvalue  of  the  one  discrete  mode, 
found  numerically  from  Eq.  16,  corresponds  to  k^= 0.125063. 

Next,  f(k)--Eq.  43— was  calculated  numerically  at  some  30  values 
of  k,  widely  spaced  where  f varied  slowly  but  more  closely  spaced  in 
the  two  dominant  resonances.  Figure  1 sketches  f(k),  from  k=0.04 
to  k^.  For  k<0.04,  f(k)  merely  falls  slowly  and  monotonically . Five 
points  at  which  cos^^H)  vanishes  were  found.  In  Fig.  1 they  are 
marked  with  circles  enclosing  the  pertinent  values  of  2n+l. 


Obviously  only  the  last  two  resonances,  appearing  toward  k=kfe, 

are  significant.  These  parts  of  the  curve  were  calculated  at  smaller 

intervals  and  the  two  maxima  were  located  with  moderate  care 

(at  kj.=0. 11482  and  4^=0.121934) . Thereafter  f(h)  was  regarded  as 

basically  a function  f (k)  to  which  the  resonances  were  added,  f 

o'  ' o 

coincides  with  f(k)  for  k<0.105,  bu<.  then  it  is  considered  to  vary 
so  smoothly  before  vanishing  at  k^  that  its  contribution  to  <PC  (Eq.  42) 
is  negligible.  Each  resonance  was  numerically  fitted  by  a gaussian, 


pj~^[f(k)  ‘ <fo>)  = a exP(-b*2) 


K is  k-kj  -j-j,  as  pertinent,  and  <fQ>  is  an  averrge  value  of  fQ(k)  over 
the  base  of  each  spike.  The  constants  were  aj=3 .551,  bj=l, 632x10^, 
aII=8.040,  and  bjj=6. 56x10^ . These  gaussians  fitted  the  resonances 
fairly  well  except  that  they  somewhat  overemphasized  the  discreteness 
of  the  "spikes"  (see  Fig.  l)  by  lowering  and  sharpening  the 
intermediate  minimum.^ 


FIGURE  1 

f(k)  versus  k in  (Eq.  43)  AT  100  Hx 
CIRCLED  NUMBERS:  ?n  ♦!  IN  * (2n  ♦l)»/2 


In  this  range  of  k,  even  for  r being  much  smaller  than  any  values 

shown  in  Stickler's  Fig.  11,  H^(kr)  is  given  with  great  precision 

by  its  asymptotic  form,  Eq.  B9>  with  x=(kT  tt+k) r.  Moreover,  (kr)1/2 

1/2 

is  very  well  approximated  by  (k^  ^r)  ' because  the  resonances  are 
narrow.  Then  the  contribution  Aj.  to  <PC  from  resonance  I is 


The  same  expression,  with  subscripts  II,  holds  for  the  other  resonance. 
The  integral  is  a standard  form  and  the  numerical  results  are  (with  r 
in  feet) 

Aj.  a 0.03669  i'1/2  r'1/2  exp(-1.532  X 10-6r2)  exp(-0.1l482ir)  ; (47) 

Aj-j  * 0.04020  i-1/2  r-1/2  exp(-3.8ll  X 10-7r2)  exp(-0.121934ir)  . (48) 

The  contribution  of  the  one  discrete  mode  is 

AM  a O.Olk&l'1/2  r"1/2  exp(-0.125068ir)  . (49) 

Equations  47  through  49  show  that  the  virtual  modes  attributable  to 

resonances  I and  II  behave  like  discrete  modes  except  for  the  real 
2 2 

exponentials  in  r . The  coefficient  of  r is  inversely  proportional 
to  the  constant  b in  Eq.  45.  Therefore,  regardless  of  a resonance's 
height,  the  virtual  mode  derived  from  a narrower  resonance  falls  off 
more  slowly  with  increasing  range.  A virtual  mode  must  decrease  more 
rapidly  with  range  than  does  a discrete  mode.  The  reason  is  seen  in 
the  integral  of  Eq.  46.  The  subset  of  continuous  modes  that  make  up 
a resonance  embraces  a small  spread  of  k and  therefore  of  k.  Regardless 
of  range,  this  spread  is  governed  by  the  gaussian  in  Eqs.  45  and  46. 

For  sufficiently  small  r,  exp(ixr)  in  Eq.  46  is  nearly  unity;  the 
modes  in  this  subset  are  substantially  in  phase.  As  r increases. 
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however,  exp(-iKr)  oscillates  in  phase  over  the  practical  spread 

of  K,  and  the  continuous  modes  interfere  even  more  destructively. 

In  the  limit  when  a resonance  is  well  approximated  by  a delta  function 

in  K,  the  consequent  virtual  mode  falls  off  only  as  just  as  a 

discrete  mode  does.  (With  their  different  fit^  to  a resonance, 

7 

Tindle  and  Guthrie  found  an  excess  falloff  in  the  form  of  a decaying 
exponential  in  r,  not  r ; however,  the  physical  cause  is  the  same.) 


A few  numerical  trials  showed  that  for  rSlOOO  ft,  all  three  modes 
interact  appreciably  as  their  phases  change  with  r.  Before  r=2000  ft 
is  reached,  only  II  and  M interact  significantly;  I has  fallen  too 
much  in  amplitude.  Eventually,  mode  II  fades  out  also,  and  a plot 
of  transmission  loss  becomes  essentially  smooth  and  monotonic. 


Numerical  values  of 
by  using  the  mode  sum 


were  calculated  at  various  values  of  r 


Three  extrema  in  the  interval  1000<r£;  5000  ft  were  located  near  the 
corresponding  ranges  in  Fig.  11  of  Ref.  3«  Long  before  r=3000  ft  is 
reached,  Aj.  has  become  negligible  |Aj-jl  is  so  small  that  the  factor 
r-^  from  |^|2  hides  any  extremum.  However,  a weak  maximum  of  r|A|2 
lay  near  r=3000  ft.  Two  still  weaker  extrema  were  found  at  greater 
ranges  by  the  differentiation  of  an  analytical  approximation  to 
r|A|2. 

These  six  ranges  are  listed  in  line  (w)  of  Table  I;  line  (S)  gives 

14 

Stickler's  corresponding  ranges. 
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TABLE  I 


COMPARISON  OF  RANGES  r (IN  feet)  AT  WHICH  EXTREMA 
OF  THE  MODE  SUM  OCCUR 


Extremum 

Min. 

Max. 

Min. 

Max. 

Min. 

Max. 

(W)  r 

1170 

1700 

2950 

3770 

4740 

5710 

(S)  r 

1120 

1740 

3000 

3820 

5420 

5660 

Discrepancy 

4.5  $ 

2$ 

2$ 

1$ 

14$ 

1$ 

Agreement  for  the  first  four  entries  is  good.  The  average  magnitude 
of  fractional  discrepancy,  2.4$,  is  quite  acceptable  in  view  of  our 
numerical  approximations  and  the  uncertainty  in  reading  ranges  from 
Stickler's  graph.  Our  last  two  extrema  are  so  weak  that  even  minor 
changes  in  the  underlying  parameters  could  alter  their  locations 
considerably.  Hence,  neither  the  good  nor  the  bad  case  is  particularly 
significant.  The  unweighted  average  value  for  magnitudes  of  fractional 
discrepancies  for  all  six  cases  is  only  4$.  The  ranges  in  line  (w) 
are  fixed  primarily  by  the  values  of  k^,  k^,  and  k^,  but  the  real 


exponentials  in  r (Eqs.  47  and  48)  have  an  appreciable  influence. 


A comparison  of  our  results  with  Stickler’s  appears  in  Fig.  2 
where  the  calculated  transmission  loss  TL  in  decibels. 


TL  = 10  log10(  H2) 


(51) 


is  plotted  versus  range  r in  feet.  The  dotted  line  corresponds  to 


const./r,  i.e.,  to  "ARL  Discrete"  of  Fig.  11,  Ref.  3,  or  simply  to 

for  +,*e  discrete  mode.  This  curve  was  used  to  make  the  vertical 
scaling  agree  with  Stickler's  graph.  The  broken  dashed  line  is  "ARL 


Discrete  Plus  Continuous"  of  Ref.  3*  The  full  line  illustrates  our 
results;  vertical  arrows  mark  the  last  two  extrema  of  Table  I.  Beyond 
r=4000  ft,  the  full  and  dotted  lines  are  essentially  indistinguishable, 
as  are  the  dashed  and  dotted  lines  for  r^7500  ft,  because  BLI  has 
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become  negligible.  Between  r=2000  and  r=2500  ft,  the  full  and  the 
dashed  curves  practically  coincide. 

The  full  line  represents  chiefly  the  sum,  with  interaction,  of 
modes  M and  II.  In  the  range  500~ri,1500  ft,  however,  the  weaker 
virtual  mode  1 contributes  appreciably.  Omitting  it  would  have 
markedly  weakened  the  agreement  with  Stickler's  result.  For  r£2500 
ft,  mode  I is  negligible  and  mode  II  has  decayed  with  range  much  too 
fast  to  yield  anything  close  to  the  dashed  line.  Replacing  Eq.  45 

13 

by  a better  analytical  fit  with  its  slower  decay  in  range  might 
improve  agreement  beyond  2500  ft. 

The  degree  of  success  out  to  rs;2500  ft,  before  the  excessive 
decay  rate  of  mode  I becomes  important,  is  gratifying  when  it  is 
remembered  that  we  have  used  only  a Pekeris  approximation  to  Stickler's 
Fig.  8.  In  particular,  the  arguments  in  Sections  V and  VI  that  only 
the  strong  resonances  need  be  included  in  evaluating  BLI  are  supported 
by  Fig.  2.  Although  no  direct  analytical  connection  has  been  shown 
with  Bartberger's  improper  modes,  we  have  at  least  found  that 
resonances  in  BLI  can  produce  similar  results.  A little  more  informa- 
tion emerges  in  Section  VIII. 

VIII.  Application  to  Stickler’s  Figure  10 

Figures  10  and  11  of  Ref.  5 differ  only  in  respect  to  frequency — 
50  Hz  for  Fig.  10  and  100  Hz  for  Fig.  11.  Stickler  found  that  a single 
discrete  mode,  near  cutoff,  still  exists  at  50  Hz.  The  two  plots  in 
Fig.  10 — "ARL  Discrete”  (AD)  and  "ARL  Discrete  Plus  Continuous"  (ADC) 
are  quite  unlike  their  counterparts  in  Fig.  11.  First,  out  to  the 
extreme  range  of  calculations  (rR90,C00  ft)  ADC  of  Fig.  10  lies  several 
decibels  above  AD  instead  of  converging  to  it  for  r-^7000  ft.  Second, 
ADC  shows  no  indication  of  mode  interaction.  Third,  Bartberger’s 
calculated  points  now  fall  on  AD,  not  ADC. 
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Numerical  trials  suggested  that  ADC  could  be  fitted  by  the  sum 
of  two  terms,  one  varying  as  (presumably  the  discrete  mode) 

and  the  other  (of  much  larger  amplitude)  varying  approximately  as 
r_1  (presumably  from  BLl) . When  a discrete  mode  is  quite  near 
cutoff,  a BLI  contribution  in  r”1  is  to  be  expected. ^A5>l6 


To  make  quantitative  comparisons,  we  proceed  as  before  by  fitting 

a Pekeris  approximation  to  Stickler's  Fig.  8.  A trial  value  of  c^, 

4985  ft/sec,  was  chosen  near  the  average  of  the  varying  c in  Fig.  8, 

and  the  eigenvalue  k,  0?  the  discrete  mode  was  calculated  from  Eq.  l6. 

The  modes  being  near  cutoff  implies  that  its  properties  are  highly 

2 

sensitive  to  the  value  of  k^  . A better  eigenvalue  was  sought  by 
using  Stickler's  actual  Cy(z)  as  the  basis  for  a first  order 
perturbation  correction.1^  The  result,  kj=0. 0621123,  differed  only 
a little  from  the  first  try.  A new  constant  value,  0^=4985. 73  ft/sec, 
was  found  to  agree  with  the  revised  k^.  Next,  the  normalization 
constant  N^  was  calculated  frcm  Eq.  l8.  Near  cutoff,  the  eigenfunction 
retains  appreciable  amplitude  deep  into  the  bottom.  Consequently,  N^ 
is  quite  small;  the  discrete  mode  amplitude  is  only  about  5$  of  that 
in  Eq.  49. 


Numerical  calculations  of  f(k)  revealed  a pronounced  spike  with  a 
peak  very  near  k=kfa.  (This  spike,  a "pseudoresonance, " is  less  than 
half  of  a true  resonance.  It  rises  smoothly  as  k-»k^  and  then  is 
abruptly  reduced  to  zero  by  the  factor  fig  in  the  numerator  of  Eq.  43. 
The  actual  resonance  condition  on  could  be  met  only  for  k>k^. ) No 
significantly  strong  resonance  was  found  at  smaller  k.  Across  the 
spike,  k varies  only  slightly,  whereas  pg  is  much  smaller  and  varies 
rapidly.  Hence,  it  seemed  best  to  change  variable  from  k to  <=(32,  as 
in  Appendix  A and  Ref.  7,  replacing  f(k)dk  in  BLI  by  g(k)<3x.  Expanding 
f(k)  in  powers  of  k showed  that  over  the  width  of  the  spike,  g is  well 
represented  by 


HESSE 


& 


s(«)  = 


aX 


(1+bK2) 


(52) 


The  constants  a and  b can  be  found  analytically,  but  as  a measure  of 
correction  for  higher  terms  emitted  from  Eq.  52  they  vere  instead 
determined  empirically  to  fit  numerical  values  of  g.  To  obtain 
known  functions  from  integration  of  g over  the  spike,  Eq.  52  was 
curve  fitted  to  another  form. 


s(<)  = (§)  t1  " expf-A*)] 


(53) 


Agreement  with  Eq.  52  was  good  (within  an  average  2 $ wherever  g 
was  of  substantial  size).  The  Hankel  function  in  BLI  was  again 
replaced  by  its  asymptotic  form  for  large  kr  (Eq.  B9),  and  k was 
expressed  (to  better  than  a part  in  10^  over  pertinent  k's)  by 


k = 


kb- 


2kb 


(54) 


Equation  52,  confined  to  the  spike,  becomes 


i - f eyiW)  A 


&K 


(55) 


Here  C*  is  a combination  of  constants  already  mentioned.  The  upper  limit 
of  k is  safely  extended  to  infinity  because  the  spike  is  narrow.  The 
expressions  in  { ) can  be  converted  to  standard  forms  in  terms  of  Fresnel 
integrals.*  The  discrepancy  £g (k)  between  Eq.  52  and  53  is  small  in 
amplitude  and  is  very  narrow.  A rough  calculation,  based  upon  a 
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triangular  fit  to  Zg,  showed  that  its  contribution  to  <p  varies 

-l/2  c 

nearly  as  r ' and  in  amplitude  is  less  than  1$  of  the  discrete 

mode. 

The  total  contribution  A to  cp  is 


A = 7.0  x ICf^ir)-1/2  expt-ikjT) 

+ 0.227  r'1  expC-UyHl  + (l+i)(f+ig)]  . (56) 


The  first  tem  comes  from  the  discrete  mode  and  the  second,  from  the 

pseudoresonance j f and  g,  functions  of  (350/r^/2),  are  defined  and 

tabulated  in  Ref.  18.  The  second  term  in  Eq.  55  is  larger  than  the 

extreme  asymptotic  result  found  through  a different  method  by 
15 

Brekhovskikh.  k^  and  k^  in  Eq.  55  are  so  nearly  equal  that  even 
for  r=9°>000  ft  the  two  exponentials  differ  negligibly  in  phase — 
another  reason  for  the  absence  of  appreciable  mode  interaction  in 
Stickler's  Fig.  10.  The  magnitude  of  Eq.  56,  converted  to  decibels, 
is  our  approximation  to  ADC  in  Fig.  10  of  Ref.  3-  The  first  term 
in  Eq.  56,  similarly  converted,  corresponds  to  AD. 

When  (as  in  Section  VIl)  it  is  assumed  that  our  AD  and  Stickler's 
are  identical,  our  calculated  ADC  differs  from  his  by  no  more  than 
0.1  to  0.2  dB,  which  is  within  the  uncertainty  of  reading  his  graph. 

It  appears  that  Fig.  10  is  adequately  explained  as  the  combination  of 
a single  discrete  mode  and  a virtual  mode  from  a pseudoresonance 
induced  by  the  discrete  mode's  being  very  near  cutoff. 

The  assumption  that  our  AD  and  Stickler's  can  be  equated  may 
seem  gratuitous  because  we  have  dealt  only  with  a Pekeris  approximation 
to  Stickler's  Fig.  8.  With  respect  to  our  Section  VII  and  Fig.  2, 
there  could  well  be  doubt.  At  short  ranges  the  one  discrete  mode  and 
the  two  virtual  ones  axe  of  comparable  amplitudes,  and  appreciable 


I 

i 

j 

i 

t 


29 


changes  in  any  one  amplitude  would  not  greatly  alter  Fig.  2.  The 
good  results  in  this  present  section,  however,  seem  to  validate  the 
assumption.  Suppose,  for  example,  that  our  Pekeris  approximation 
drove  the  discrete  mode  too  close  to  cutoff.  Then  g^H  is  pushed 
downward  too  close  to  the  value  n/2.  Our  pseudoresonance  and  the 
virtual  mode  are  made  too  strong.  Simultaneously,  the  discrete  mods 
is  made  too  weak  because  the  term  |tan(gjH)|  in  Eq.  18,  which  strongly 
dominates  the  value  of  becomes  too  large  and  makes  too  small. 

, Therefore,  our  calculated  differences  between  AD  and  ADC  would  have 

been  decidedly  larger  than  Stickler's  differences.  A similar 
argument  would  hold,  of  course,  if  our  discrete  mode  had  been  too 
far  from  cutoff;  our  AD  and  ADC  would  have  been  too  close  together. 

Two  speculative  comments  might  be  added.  The  first  concerns  the 
abrupt  snift  of  Bartberger's  calculated  points  from  lying  on  AD  in 
Fig.  10  of  Ref.  3 to  lying  on  ADC  in  Fig.  11.  It  would  appear  that 
for  some  reason  not  yet  clear,  the  method  of  improper  modes  can 
duplicate  the  effects  of  true  resonances  within  BLI  and  the  consequent 
virtual  modes,  but  that  it  does  not  detect  a pseudoresonance  caused  by 
a discrete  mode  very  near  cutoff.  (The  method  would,  of  course,  recognize 

| 

and  use  that  discrete  mode  as  well  as  any  others  that  might  be  present.) 

The  second  comment  concerns  Kutschale's  computed  points,  shown 
in  Fig.  10  of  Ref.  3.  These  points  might  be  interpreted  as  showing 
mode  interaction  between  the  discrete  moue  and  BT-T — a small  effect  in 
amplitude,  but  much  larger  than  any  found  by  Stickler  and  by  us.  One 
implication  would  then  be  that  Kutschale's  discrete  mode  and  BLI  are 
much  nearer  in  amplitude  than  we  have  found . A more  definitive  test 
is  afforded  by  the  corresponding  interaction  distance  Ar  (defined  by 
Ar=&t/Ak,  where  in  our  calculation  Ak  is  k^-k^).  Kutschale's  Ar  can 
be  roughly  estimated  as  45,000  to  60,000  ft;  in  sharp  contrast,  our 
Ar  would  be  about  8x10  ft.  This  seeming  interaction  may  be  only  an 
accidental  result,  but  in  any  event,  it  shows  that  different  computer 
methods  can  lead  to  noticeably  different  results . 
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IX.  Summary 


Under  the  restrictions  of  constant  water  depth,  only  vertical 
variations  in  parameters,  no  absorption,  and  a bottom  that  can  be 
treated  as  a fluid  half-space,  Sections  I through  IV  give  the  nr>wn»i 
mode  method  for  propagation  of  underwater  sound.  The  language  of 
eigenfunctions  and  eigenvalues  is  used,  with  no  recourse  to 
integration  in  the  complex  plane.  Only  an  omnidirectional  point 
source,  cw  at  a single  frequency,  is  used.  Except  perhaps  for  the 
two  postulates  in  Section  I,  nothing  new  is  presented;  the  goal  was  a 
compendium.  For  the  Pekeris  problem,  it  was  demonstrated  that  the 
integral  over  the  “continuous”  modes  is  identical  to  the  branch  line 

p 

integral  (BLl)  found  by  Ewing,  Jardetzky  and  Press  through  use  of 
integration  in  the  complex  plane.  (Tindle  and  Guthrie  have 
independently  and  simultaneously  derived  the  same  result.)^  Since  the 
Btl  of  Bef.  2 differs  from  that  of  Pekeris,1  the  outcome  is  of  interest. 

Sections  V and  VI,  supplemented  by  Appendices  A and  B,  discuss 

properties  of  the  Btl,  i.e.,  the  integral  over  continuous  modes,  for 

the  Pekeris  problem.  It  is  proved  analytically  that,  piecewise,  BLI 

contributes  practically  negligibly  to  the  sound  field  except  where  the 

integrand  displays  narrow  "spikes"  or  resonances  (see  Fig.  l).  The 

conditions  for  such  resonances  to  exist  with  the  consequent  production 
6 7 

of  virtual  modes  ’ ' are  demonstrated.  (Again,  Ref.  7 covers  much  of 
the  same  ground  although  with  a different  but  equivalent  integral  and 
with  less  analytical  proof.) 

Finally,  in  Sections  VII  and  VIII  the  preceding  methods  are  used 
to  afford  comparisons  with  some  computer  calculations  by  Stickler^ 
(specifically,  his  Figs.  8,  10,  and  11 ).  For  this  purpose,  his  assumed 
model  was  approximated  by  a Pekeris  problem.  Reasonably  good  agreement 
is  shown  with  his  Fig.  11  (our  Table  I and  Fig.  2)  and  excellent 
agreement,  with  his  Fig.  10.  Advantages  of  the  semianalytic  method 
used  here  are  that  (l)  a partial  separation  of  range  variation  is 
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effected  and  (2)  fairly  simple  explanations  are  afforded  of  su'ii 

considerable  differences  as  those  between  Stickler's  Fig.  10  and  ! 

Fig.  11.  | 

i 

t 

% 

An  empirical  conclusion  about  Bartberger's  method  of  improper  ; 

4 5 | 

modes  * is  presented;  it  is  based  only  upon  specific  results  found  j 

here.  That  method  seems  able  to  duplicate  the  effect  (virtual  modes) 

of  true  resonances  in  BLI  but  not  to  handle  the  behavior  of  BLI  when 

a discrete  mode  is  near  cutoff. 
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At  the  beginning  of  Section  V it  was  suggested  that  the  integrands 
of  Eqs.  32  and  33  might  have  singularities  at  k=k^,  where  Pg  vanishes, 

provided  also  that  S1H=(n+l/2)«  there.  We  now  show  that  no  singularity 

J-  2 2 

exists.  Using  Eqs.  14-  and  19  for  P^  and  pg  (with  k^  replaced  by  k ) 

and  abbreviating  k^-k^2  £is  K2  and  Pg  as  k,  we  have 

0 o l/2 

px  = + (1 r+»r)  . (Al) 


2 2 

In  Eq.  32  we  multiply  numerator  and  denominator  by  6 =( p^/p^)  and  at 
first  consider  only  the  resulting  denominator  D(k),  assuming,  of  course, 
that  neither  sin(P1d)  nor  sin(p-jZ)  vanishes  at  k=kb: 

D(k)  = 52Px2  - (62  ^-K2)  sin^H) 

= 62Px2  - (&2P12-/C2)  cos2©  ; (A2) 

0 = ( n+l/2 )n  - PXH  . (A3) 

These  equations  are  exect . We  now  assume  that  k is  so  near  k^, 
where  both  K and  cos(P1H)  vanish,  that  k^«K*  and  9^«(n/2)^ . The 
consequences  are  that 


cos2©  = 1 - 02  + . . . ; PXH  = KH  + (l^)^2  + • • • • (A4) 

The  first  omitted  terms  are  in  ©^  and  k*  . Since  P^H-> (n+l/2  )it  as  K-*0, 
evidently  KH  equals  (u+l/2)n.  Equations  A3  and  A4  yield 


© = - 


t 0(Kk) 


(A5) 
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Finally,  Eqs . A2,  Ak,  and  A5  show  that 


■i 

b; 

i 


D(k)  = 2[ 1 + 0(k2)J  . (A6) 


Returning  to  Eq.  32,  we  see  in  the  numerator  a factor  ^=K  anc*  the 
quantity  kdk.  Since  fi^=K2  is  k^-k2,  it  follows  that  kdk=-/cdx. 
Accordingly,  2 cancels  in  numerator  and  denominator.  The  integrand 
of  Eq.  32  does  not  vanish  at  k=k^,  but  remains  finite. ■ 

Equation  33,  which  we  have  not  used,  behaves  itself  similarly. 

In  comparison  to  Eq.  32,  it  lacks  the  factor  fL  in  the  numerator,  but 

1/2  o 

its  denominator  becomes  D ' « K,  instead  of  D . 
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a rt 


. 


is  known^: 


The  pertinent  integral  for  ^ 

J,  -(r)l1  ho2)W  * r [HoWhi2)W 


(Bl) 


- Hi(x)Hi2)w 


fr 


x=ki-ir 


HB2 

's  are  Struve  functions.  Continuing  with  the  abbreviations 

x=kr,  x^=k^r,  and  x^  ^k^  r,  we  write 


H0(x)  = Y0(z)  + P(x)  5 Hi  “ \(x)  + «(x)  • (®2) 


Equations  B2  are  unexceptionable  provided  that  P and  Q are  defined 


properly;  we  choose  this  form  with  an  eye  on  later  use  of  asymptotic 


expressions.^  Substituting  Eqs . B2  in  Eq.  Bl  and  using  the  definitions 


of  Hankel  functions. 


t(2). 


H'  {(x>  = J -j (x)  - iY  _(x) 

0,1'  ' 0,1'  ' 0,1'  ' 


(B3) 


we  obtain 


J , ■ (¥)  Ho2)<*>  - 1 tJ<,WYiW  - Jl^)Yo(x» 

+ P(x)h[2)(x)  - (|)  Q(x)H^2\x) 


(B4) 


‘1-1 


B4 


The  [ ] in  Eq.  B4  has  the  value  (2/jtx).  Let  us  express  P and  Q in 
terms  of  p(x)  and  q(x): 


41 


2Sn 


5 


P(x)  = (“^)  u + p(x)]  ; Q(x)  = [1  + q(x)]  . (B5) 


Equation  B4  reduces  to  the  form 

J? I = (^r)  [1  + H3_2)(x)  + p(x)H^2^(x)  - xq(x)H^(x)]  I 1 • (b6) 

xi-i 

it 

t When  is  evaluated  at  its  limits  and  the  result  is  f added  to 
Jt  (Eq.  4l)f  and  after  Eq.  37  is  used,  we  find  for  JMr 

J) t = r'1  [f^H^^r)  - f (lt£_1)HJ2 ^ _1r ) ] 


[p(kir)H[2)(kJr)  " P(k/_1r)Hl2)(k|-lr)} 


(B7) 


tkfr  q(kir)H^2)(k/r)  - k^r 


(2)/ 


The  next  step  is  to  sum^ from  Z-1=0  (at  k=0)  to  corresponding 
to  ki^k^.  The  first  bracket  in  Eq.  B7  leads  to  a telescoping  sum 
that  reduces  to  its  first  term  (k=0)  and  to  its  last  (k=kb).  We 
recall  from  the  comment  after  Eq.  41  that  a term  -2i/nr  is  to  be 
ignored.  But  f(o)  vanishes,  and  so  does  f(k^)  except  in  the  unlikely 
event  discussed  at  the  start  of  Section  V.  We  disregard  this  special 
case  (see  Appendix  A).  Accordingly,  the  first  bracket  in  Eq.  B7  sums 
to  zero. 


£j 

The  remaining  two  terms  of  Eq.  B7  in  { ) are  nominally  to  be  j 

summed  starting  with  1-1=0,  i.e.,  k=0  and  1=1.  However,  in  the  first  \ 

interval  k to  k,  the  function  f(k)  vanishes  at  k=0:  therefore  a,  is  i 

o 1 v ' 1 

identically  zero.  The  sum  is,  then,  from  1=2  to  A.  It  appears  safe 

® B5 

to  replace  p(x)  and  q(y)  by  their  first  asymptotic  terms  J (respectively 
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2 2 2 
-l/x  and  +l/x  ),  which  fall  off  as  l/r  . The  next  asymptotic  terms 

. L. 

go  much  faster,  as  l/r  . In  contrast,  each  discrete  mode  falls  off 
asymptotically  only  as  r~‘L/^.  The  summand  in  Eq.  B7  is  now 


Jl \ = -(—)  (yr^^v) + (v)~M2)(k*r) 

+(r)  (ki.ir)'2Hi2)(ki.ir) + (ki-ir)“M2)(Vir) 


(B8) 


(2) 

At  this  point  we  replace  the  Hankel  functions  IP  ' by  the  first 
terms  of  their  asymptotic  expressions: 


■ «> 

For  the  same  arguments  x,  H^2^  and  have  equal  magnitudes,  but  the 
terms  in  H£ ' are  multiplied  by  l/r^,  not  l/r  as  for  those  in  IP  1 . 
Hence,  we  retain  only  the  latter  terms.  (Discarding  H^2^  is  not 
essential.  It  could  be  retained  and  handled  as  we  proceed  below.  In 
that  event,  however,  the  first  term  omitted  from  Eq.  B9  should  also 
be  used.)  Equation  B8  becomes 


eW"*  r-5A  a /k-3/2 


(BIO) 


When  the  portion  of  Sq.  BIO  that  depends  upon  l is  summed,  it  is  found 
to  telescope  in  part,  leaving 


( J 


-ik.  r 

-a#k-V2e  ^ + 
b b 


-1 


Yj  <Vl-al>ki 


-5/2 


-ik^r 


. (Bll) 


StiVi! 


From  Eq.  38  it  is  evident  that 


- a*  * ' 


1 £ l £ L 


(BI2) 


Fran  the  same  equation,  it  follows  that  f(k^)=m^  k^+a^  at  f=i^, 
but  there,  ordinarily,  f vanishes.  Hence  we  can  replace  a^  by 

-m*bkV 

Restoring  the  factor  -2i( p^/p^)  omitted  after  Eq.  37  and.  recalling 
that  we  have  calculated  the  negative  of  BLI,  we  obtain  finally 


„ = 215/2(f»V2\l/2  r-5A 

c Vb /W 

b 1 

m m 

which  is  Eq.  42. 

Note;  Bart  of  the  analysis  here  and  in  Section  V depends  on  the 
assumption  that  f(k)=M|k|  at  sufficiently  small  real  or  imaginary  k. 
Near  k=0,  0^  is  nearly  equal  to  k^  and  0g,  to  k^.  The  denominator 
of  f(k)  (Eq.  43)  can  cause  no  trouble  even  if  cos^^)  or  sin(p^H) 
should  vanish;  condition  2 of  Section  VI  is  not  met.  However, 
sin(P1d),  sin(0.jz),  or  both  might  vanish  at  k=0.  We  outline  the 
analysis  for  vanishing  sin^d). 

p O 

With  k «ky  , whether  k is  real  or  imaginary,  Eq.  14  yields 


(B13) 


id  = k d - 
K1  w 


I3T7 


(B14) 
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